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Abstract 

We estimate the sizeable cross section for deep exclusive electroproduction of an exotic J PC = 
1 h hybrid meson in the Bjorken regime. The production amplitude scales like the one for 
usual meson electroproduction, i.e. as 1/Q 2 . This is due to the non- vanishing leading twist 
distribution amplitude for the hybrid meson, which may be normalized thanks to its relation 
to the energy momentum tensor and to the QCD sum rules technique. The hard amplitude 
is considered up to next-to-leading order in as and we explore the consequences of fixing 
the renormalization scale ambiguity through the BLM procedure. We study the particular 
case where the hybrid meson decays through a irr] meson pair. We discuss the ttij generalized 
distribution amplitude and then calculate the production amplitude for this process. We 
propose a forward-backward asymmetry in the production of tt and 77 mesons as a signal for 
the hybrid meson production. We briefly comment on hybrid electroproduction at very high 
energy, in the diffractive limit where a QCD Odderon exchange mechanism should dominate. 
The conclusion of our study is that hard electroproduction is a promissing way to study exotic 
hybrid mesons, in particular at JLAB, HERA (HERMES) or CERN (Compass). 
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1 Introduction 



Within quantum chromodynamics, hadrons are described in terms of quarks, anti-quarks and gluons. 
The usual, well-known, mesons are supposed to contain quarks and anti-quarks as valence 3 degrees of 
freedom while gluons play the role of carrier of interaction, i.e. they remain hidden in a background. On 
the other hand, QCD does not prohibit the existence of the explicit gluonic degree of freedom in the form 
of a vibrating flux tube, for instance. The states where the qqg and gg configurations are dominating, 
hybrids and glueballs, are of fundamental importance to understand the dynamics of quark confinement 
and the nonperturbative sector of quantum chromodynamics PQ-[I]. 

The study of these hadrons outside the constituent quark models, namely exotic hybrids, is the main 
reason of the present paper. We investigate how hybrid mesons with J PC = 1 ^ may be studied through 
the so-called hard reactions. We focus on deep exclusive meson electroproduction (see, for instance 
0) which is well described in the framework of the collinear approximation where generalized parton 
distributions (GPDs) [x] and distribution amplitudes 8 describe the nonperturbative parts of a factorized 
amplitude HJ. 

In a previous paper we showed that contrarily to naive expectations, the amplitude for the 
electroproduction of an isotriplet exotic meson with J PC — 1 h may be written in a very similar way as 
the amplitude for non-exotic vector meson electroproduction. The main observation of our work was that 
the quark-antiquark correlator on the light cone includes a gluonic component due to gauge invariance and 
leads to a leading twist hybrid light-cone distribution amplitude. In this paper, we extend our analysis 
of the electroproduction process and calculate the differential cross section as a function of Q 2 . 
We also study the hybrid meson as a resonance in the reaction ep -> ep(ir°rj). The first experimental 
investigation of the hybrid with J PC = 1 h as the resonance in 7r — 77 mode was implemented by the 
Brookhaven collaboration E852 [^. Present candidates for the hybrid states with J PC = 1 h include 
7Ti(1400) which is mostly seen through its irq decay and 7Ti(1600) which is seen through its nrj' and irp 
decays ^2- Theoretically these states are objects of intense studies £Q, mostly through lattice simulations 




2 Hybrid meson production amplitude 

We propose to study the exotic hybrid meson by means of its deep exclusive electroproduction, i.e. 

e(fei) + N(jn) -» e(fc 2 ) +H{p) + N(p 2 ), (1) 

where we will concentrate on the subprocess: 

ll{q) + N{p x ) -> H L (p) + N( P2 ) (2) 

when the baryon is scattered at small angle. This process is a hard exclusive reaction due to the transferred 
momentum Q 2 is large ( Bjorken regime). Within this regime, a factorization theorem is valid, at 
the leading twist level, which claims that a partonic subprocess part described in perturbative QCD 
(pQCD) can be detached from universal soft parts, which are generalized parton distributions and meson 
distribution amplitudes. Below we will analyze in more details how this factorization theorem applies to 
the process under study. 

In this paper, the main object of our investigation is the isotriplet of mesons with quantum numbers 
J PC = 1 '". Such mesons can be named as exotic mesons because they do not exist within the usual 

3 The valence degrees of freedom define the charge and other quantum numbers of corresponding hadrons, while the sea 
configurations do not change the quantum numbers. 
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quark model. To illustrate the latter we shortly remind the main steps of the description and classification 
of meson states in the quark model. 



2.1 Quark model and spectroscopy 



It is well known that in the quark model the hadrons, mesons and baryons, are bound states of quark- 
antiquark or three-quarks systems. Let us consider the mesons, i.e. the quark-antiquark systems. Their 
total angular momentum results from the summation of spin 5 and orbital L angular momenta of quarks. 
Neglecting a spin-orbital interaction, the quantum numbers 5 and L may be considered as additional 
quantum numbers for the classification of hadron states. Therefore, the eigenvalues of the squares of the 
angular momenta read: 

J 2 = J(J+1) S 2 = 5(5 + 1) L 2 =L(L + 1), 

J = S + L, (3) 

where the number L may take all positive integer values (including zero). The meson octets correspond 
to the case where S = Q, 1. For given values of 5 and L, the total angular momentum J can take the 
values 



J = S + L, S + L-l, ...,\S-L\. (4) 

The values 5 and L are related to the C- and P-parity of the quark-antiquark system in the form: 

C = (-) L+S , P=(-) L+1 . (5) 

Consequently, in the quark model, the quantum numbers 5, L, J, P, C and the relations between them 
(0J , (J2J) allow to introduce the following classification of the meson states : 

• 5 = 0, L = J : 

J PC = Cr+, l+-,2-+ 3+-, ... (6) 



S = 1, L = 0, J= 1 



J 



PC 



(7) 



1, L = 1, J 



2, 1,0 



J PC = 0++, 1 



-+ o++ 



(8) 



S = 1, L = 2, J = 3, 2, 1 



J PC = 1- 



, 2—, 3" 



(9) 



and so on. From this, one can see the mesons with J PC = 0—, 0+-, 1"+, ... , are forbidden. However, 
such mesons may be described beyond the quark model. Indeed, we may add an extra degree of freedom 
(a gluon, for instance) to get the needed quantum numbers, see for instance PJ. Below we will consider 
this case in more details. 
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Figure 1: Typical diagram describing the electroproduction of a meson at lowest order. The grey blobs are 
non-perturbative matrix elements, namely the meson distribution amplitude and the nucleon generalized 
parton distribution 



2.2 Kinematics and leading order amplitude 



Let us fix the kinematics of the deep electroproduction process. We are interested in the scaling regime 
where the virtuality of the photon Q 2 = —q 2 is large and scales with the energy of the process. We 
denote by pi (p2) the momentum of the incoming (outgoing) nucleon, while p is the momentum of 
the longitudinally polarized hybrid meson of mass Mr- We construct the average momentum P and 
transferred momentum A: 

P2+PI . a2 . nn\ 

P = - , A =;p2- pi, A z =t. (10) 

It is convenient to introduce the following two light-cone vectors: 

n* = A(1,0 T ,1), n=^-(l,0 T ,-l), n* ■ n = 1, (11) 

where A is an arbitrary dimcnsionful parameter 4 . Then the Sudakov decompositions for all the relevant 
momenta take the form: 

A M = -2^ + (M 2 n^ + Aj, A T • n — A T • n* = 0, 
-2 



M 2 2 V— A 2 

P, = n; + —n„ P=M, £<-_=^<l, 

op * 1 Q 

p,=q,-A, = 2(£ - On; + - (M^n, A£. (12) 



Here, the parameters £ and £ are related by 



M 2 H =4($-i)(^l-m Z )+A 2 T . (13) 



4 Within the infinite momentum frame, the parameter A may be chosen as P . Here, A+ = (A + A 3 )/\/2. 
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The photon longitudinal polarization vector can be written as 



Hi 

Q 



Q 

—si 

4£ 



where the notation Q = ^/ Q 2 is introduced. 

The leading order amplitude for the process corresponding to the diagrams in Fig. ^is 



A 



where the ^-matrix is given by 



dr,e^(N(p 2 )H(p)\ 



SS 



\N(pi)), 



S , = Texp|z j (Px\Cqcd{x) + £qed{x)^ 



(14) 



(15) 



(16) 



Applying the factorization theorem, this amplitude can be written at leading twist and when —t <C Q 

as 



A = J dz J dx $ H {z, p 2 , (1%) H(x, z,Q 2 , p 2 , p 2 R ) F(x, p? , Mr) = $h <8) H ® F, 



(17) 



where the parameters p, 2 and p 2 R are the factorization and renormalization scales, respectively. Through- 
out this paper, we will adopt the convention that \i = /j,r 5 . In Ea. lfT7jl . H is the hard part of amplitude 
which is controlled by perturbative QCD. The hybrid meson distribution amplitude &h describes the 
transition from the partons to the meson, and F denotes generalized parton distributions which are 
related to nonperturbative matrix elements of bilocal operators between different hadronic states. 
More precisely, the factorized amplitude (|17fl may be written as : 



where 



A 



(?) 



e^a s fHC F 
V2N C Q 



(18) 



H 



ff 



V (M, 



dx 



U(p 2 )nU( Pl )H ff ,(x) + U(p 2 ) a » a 



2M 



U( Pl )E ff ,(x) 



± 



x + £ — ie x — £ + it 



dy<t> M {y) 



1 1 

± 



y i-y 



(19) 



Here, functions H and E are standard leading twist GPD's and their properties are fairly well-known 
(see, for instance, a review of Diehl in [7j). In l|19fl . we include the definition of H~ff and y( Af >+) which 
will be useful for the comparison with the p meson case. The hybrid meson distribution amplitude which 
enters Eqn. (|19fl is a new object and we will carefully study it in the next subsection. Note that the 
simple pole over y in l|18fl does not lead to any infrared divergency if the function 4> H (y) vanishes when 
the fraction y goes to zero or unity. 



5 The arguments in favour of such a choice in the case of the pion form factor are discussed e.g. in 1131 1141 
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2.3 Hybrid meson distribution amplitude 



In this subsection, we will consider in detail the properties of the hybrid meson distribution amplitude 
(see also JO])- The Fourier transform of the hybrid meson -to-vacuum matrix element of the bilocal 
vector quark operator may be written as 



{H(p, \)$(-z/2)^[-z/2; z/2]1>(z/2)\0) = (20) 
if H M H 



' 

where with A = (0, +1, —1) describes the polarization states of the hybrid meson. It is convenient 
to define the four-vector corresponding in the ultrarelativistic limit to the longitudinal polarization 

as 

For the longitudinal polarization case, only the term with <f>^ contributes, so that 

l 

(H L (p,0m~z/2)^[-z/2;z/2Wz/2)\Q) = if B M B e<® J dye^-^ z ' 2 <t>«{y) (22) 

o 

where y = 1 — y and H denotes the isovector triplet of hybrid mesons; fn denotes a dimensionful coupling 
constant of the hybrid meson, so that (f> H is dimensionless. We will discuss its normalization later. 
In l|21l) and (|22|l . we insert the path-ordered gluonic exponential along the straight line connecting the 
initial and final points [z\ ; Z2] which provides the gauge invariance for bilocal operator and equals unity 
in a light-like (axial) gauge. For simplicity of notation we shall omit the index L from the hybrid meson 
distribution amplitude. 

Although exotic quantum numbers like J PC = 1 h are forbidden in the quark model, it does not prevent 
the leading twist correlation function from being non zero. The basis of the argument is that the non- 
locality of the quark correlator opens the possibility of getting such a hybrid state, because of dynamical 
gluonic degrees of freedom arising from the Wilson line. This may be seen easily through a Taylor 
expansion of the non-local correlator 

(H(p, 0)$(-«/2) 7/4 [-*/2; z/2}iP(z/2)\0) = (23) 

E ^■■ z »n(H(p,o)moh» B w .. B Mn v(o)|o), 

n odd 

where is the usual covariant derivative and D^= \{D^ — D^). The term corresponding to z = 
refers to the standard quark model contribution, which is zero for the exotic hybrid quantum numbers. 
The first non zero contribution arises from the first derivative contribution in this expansion, and one 
can check that, more generally, only odd terms contribute in this expansion. It is clear that due to gauge 

invariance, such occurrence of operators naturally provides gluonic degrees of freedom, which enables 
the production of hybrid state at twist 2 level. One can check explicitly that the corresponding quantum 
numbers are indeed the one of the hybrid state (see detailed discussion in JOI)- Using charge conjugation 
invariance of H°, one can show that the corresponding distribution amplitude is antisymmetric, namely 

d ) H {y)^^ H {l~y). (24) 

The result of this analysis is that the hybrid light-cone distribution amplitude is a leading twist quantity 
which should have a vanishing first moment because of the antisymmetry property of the distribution 



G 



amplitude. This distribution amplitude obeys usual non-singlet evolution equations 8; and has an asymp- 
totic limit PH] 

$ H = 30y(l-y)(l-2y) (25) 

The normalization factor (coupling constant) fjj is defined through the matrix element of the energy- 
momentum tensor |16j . It may be related, by making use of the equations of motion, to the matrix 
element of quark-gluon operator and estimated with the help of the techniques of QCD sum-rules |17|. 
One of the solutions corresponds to a resonance with mass around 1.4 GeV and normalization factor 6 

f H ~ 50MeV . (26) 

If it turns out that only one resonance can be attributed to such an hybrid state, this QCD sum-rule 
analysis is sufficient to fix the value of fjj. If the scenario with two resonances is confirmed, this value of 
fn corresponds to an effective coupling to the total contribution of these two resonances. Despite the fact 
that QCD sum rules cannot distinguish at the level of the coupling between two very close resonances 
and a single one, if one would define fx (resp. /a) the coupling to the first resonance (resp. second), one 
should write 

f H = fi + ft ■ (27) 
Thus, when selecting experimentally each resonance by their decay modes, we know for sure that one 
of the coupling should be larger than /h/v2. Thus, the fact that the exotic hybrid quantum numbers 
could be attributed to two very close resonances does not spoil the conclusion about the expected order 
of magnitude of the hybrid distribution amplitude. From now on, we will consider the case where only 
one of the 7Ti candidates is an exotic hybrid meson. 

The coupling constant fjj is the subject of evolution given by the formula, see e.g. 

fas(Q 2 )\ K ° v _2 7QQ (0) 



\a s {M z H )J f3 

where the anomalous dimension 7qq(0) = 16/9 and /3o = 11 — 2nj/3. The exponent K a is thus a small 
positive number which drives slowly to zero the coupling constant fu{Q 2 )- Since experiments are likely 
to be feasible at moderate values of Q 2 , we neglect this evolution. 



3 Cross-sections for hybrid meson electroproduction 



In this section, we focus on the computation and analysis of the differential cross section for longitudi- 
nally polarized hybrid meson electroproduction. The experimentally accessible differential cross section 
corresponds to the process Q for which the reaction © is the subprocess. We assume that the virtual 
photon has a longitudinal polarization, so that the factorization theorem is valid owing to the absence of 
infrared divergences. We restrict ourselves to the leading twist contributions. 

When we estimate hybrid meson electroproduction cross section, we systematically compare it with the 
similar contribution (i.e. without gluon GPDs) to the cross section for longitudinally polarized p meson 
electroproduction. The unpolarized cross section corresponding to the reaction is defined by 7 



d<JL 1 



dt 16n{s-m 2 N )\(s,-Q 2 ,m 2 N )2 



!EIA*)| 2 > ( 29 ) 



pol. 



where the amplitude A./q) is determined by 1)18(1 : s, t are the usual Mandelstam variables and tojv is the 
nucleon mass. The function A is standardly defined by 

A 2 (x, y, z) = x 2 +y 2 + z 2 - 2xy - 2xz - 2yz. (30) 



6 Our fjj corresponds to 2v / 2/h in the notations of Ref. 1171 . 
7 The flux factor is chosen as in I2(>) . 
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To calculate the cross section (|29ll . we need to model the corresponding GPD's. We apply the Radyushkin 
model |22] where the function H, see (TT%1) . is expressed with the help of double distributions F q (x,y;t). 
We have 



mm (ir— I 



o 

min{^,i±|} 

^ " ' ^(L^zM tVlt] . (3i) 





where standard notations are used. For the double distribution F q (X, Y;t), we assume the ansatz sug- 
gested by Radyushkin \'27\ : 

F q (X,Y;t) = j| g W6 Y( [r-X)3 y) ' < 32 > 

and a similar expression for the anti-quark contribution. As shown in [2Hj , this definition of the double 
distribution is not completely compatible with the structure of the corresponding matrix elements; intro- 
ducing D-terms restores the self-consistency of this representation. Taking into account these D-terms, 
the GPD's iJjJTJ are modified into : 

K(x,t,t) = H"(x,tt)+e^-\x\) D{x ^ t) , (33) 

where D(x/£,0) is given as in In the present paper, we concentrate on the region where the values 
of the skewedness parameter £ are rather small. Hence, it is legitimate to neglect these D-terms in the 
amplitude of p meson production. Meanwhile their contributions to the hybrid meson production always 
vanish owing to the anti-symmetric properties of these D-terms. 

In (|32l) . functions q(x) and q(x) are the ordinary quark and anti-quark distributions in the nucleon for 
which we use the MRST98 parameterization |29| . An important aspect of each model of GPD's is its 
dependence ot t j3U]. Here it is assumed to be factorizable through the functions Ff (i) for each flavour, 
which are equal to 

F 1 " = 2Ff + Ff, F 1 rf = 2F 1 n + Ff, (34) 

where Ff and Ff are the proton and neutron electromagnetic form factors. Note that we neglected the 
strange form factor because it is small. In the same way, we can write the expression for the function E. 
We neglect here its contribution because it is small and quite model-dependent. 

Finally, to get prediction for the cross sections we need to fix the renormalization scales. In order to 
estimate theoretical uncertainties of this procedure we fix the scale p, 2 R in two different ways: firstly, in 
the naive way, by assuming p R = Q 2 , and secondly, by applying the BLM prescription [2~2], 
The resulting differential cross sections for hybrid meson and p meson (quark contribution only) produc- 
tion are shown on Fig. |2]for xb — 0.18 and 0.33, using the above mentioned naive scale fixing. 
The BLM procedure, which is discussed in details in |321 ? leads to the following values of the renormal- 
ization scales: 

= e Q , tor p meson, 
H 2 R = e- 513 Q 2 , for H meson. (35) 

for the case £ = 0.2 (or ig w 0.33), and 

p R = e~ 4e8 Q 2 , for p meson, 

fi% = e~ 5 -°Q 2 , for H meson. (36) 
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Figure 2: Differential cross section for p and hybrid meson production with the naive choice of the 
rcnormalization scale and different xb- 
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Figure 3: Differential cross-section for exotic hybrid meson electroproduction (dashed line) with p? R = 
e" 513 Q 2 compared with the quark contribution to p° electroproduction (solid line) with y? R = e~ 4 9 Q 2 , 
as a function of Q 2 , for xb ~ 0.33. The dash-dotted line is the result of Vanderhaegen et al f° r P 
electroproduction. 
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for the case £ = 0.1 (or k 0.18). 

Note that, taking into account the D-terms, the p meson BLM scale is slightly diminished. For instance, 
in the case xb ~ 0.33 we have 

/4=e- 5 ' 4 Q 2 . (37) 

These renormalization scales have rather small magnitudes. This has a tendency to enlarge the cross 
sections but may endanger the validity of the perturbative approach. However, it is possible that the 
coupling constant as stays below unity and the perturbative theory does not suffer from the IR divergen- 
cies. We will use the Shirkov and Solovtsov's ansatz j25j where the analytic running coupling constant 
takes the form: 



a s (Ma) = 75- 



1 A 2 

■1 iy QCD 



]nfj,%/A 2 QCD 



(38) 



Here Aqcd is the standard scale parameter in QCD. The second term in (|38|l assures the absence of a 
ghost pole at n 2 R = Aq CD and has a nonperturbative source. Detailed discussion on this point may be 
found in |23| and references therein. 

Recently, in |2fi| the role of power corrections due to the intrinsic transverse momentum of partons (the 
kincmatical higher twist) has been investigated. In that approach the inclusion of the intrinsic transverse 
momentum dependence results in a rather strong effect on the differential cross-section before the scaling 
regime is achieved. In the renormalization scale p? R is defined by the gluon virtuality so that the 
scale is a function of parton fractions flowing into the corresponding gluon propagator. 
On Fig. |3 we present our results for the differential cross section of the hybrid meson electroproduction 
compared to the p meson electroproduction, using the BLM scales. We can see that the hybrid cross 
section is rather sizeable in comparison with the corresponding p meson cross section. We also show the 
results obtained in for the p meson electroproduction. We see that in the region Q 2 ~ 5 — 10 GeV 2 
the size of the p meson cross section obtained with the inclusion of transverse momentum effects is very 
close to the analogous cross section computed with the BLM scale and without the intrinsic transverse 
momentum dependence. On the other hand, for higher values of Q 2 the leading order amplitude computed 
with the BLM scale fixing is falling faster that the corresponding amplitude derived in Ref. > whereas 
for smaller values of Q 2 it is larger than that prediction. We do not want to claim here that kinematical 
higher twist contributions have no effects at low values of Q 2 but rather that a rather strong effect on 
the Q 2 dependence of the cross sections may be dictated by another mechanism which is much more 
controllable since it depends on the estimate of higher order perturbative contributions. 
All this shows that the scale fixing ambiguities lead to a non negligible theoretical uncertainty on the 
absolute value of cross sections. It is important however to understand that most of this uncertainty does 
not apply to ratios of cross sections, and in particular to the most interesting ratio da H : d(j p , which 
measures the expected cross section for hybrid production with respect to the well measured and large 
cross section for p meson production. Indeed, as shown on Table 1, this ratio is very insensitive to the 
scale fixing procedure. Moreover it is not small when xb is large enough and almost Q 2 independent. 
The decreasing value of the ratio when ig diminishes comes from the relative sign of the two terms 
contributing in ifTSJl. i.e. when £ — * the structure H~ goes to zero too. 

In conclusion of this section, we would like to stress that the present work has demonstrated the feasibility 
of hybrid meson production experiments in electroproduction at moderate energies. An obvious remaining 
question is how much of the cross section is observable in a dedicated experiment. If the experiment is 
able to detect the final state electron and baryon and to measure their momenta with good accuracy, 
a missing mass analysis may allow to identify and study all decay channels of the hybrid. In the next 
sections, we discuss the cases where the hybrid meson is detected through a particular decay channel. 
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X B 


0.33 


0.18 


Q 2 (GeV 2 ) 


3.0 


7.0 


11.0 


17.0 


3.0 


7.0 


11.0 


17.0 


= Q 2 


0.123 


0.123 


0.123 


0.123 


0.0325 


0.0326 


0.0326 


0.0326 




0.131 


0.133 


0.133 


0.134 


0.0356 


0.0362 


0.0365 


0.0367 



Table 1: Ratio da : da p for both the naive and BLM scales and for the different values of xb- 



4 Study of hybrid mesons via the electroproduction of tttj pairs 

In the case where there is no recoil detector which allows to identify the hybrid production events through 
a missing mass reconstruction, one will have to base an identification process through the possible decay 
products of the hybrid meson H°. Since the particle 7Ti(1400) has a dominant nr) decay mode, we now 
proceed to the description 8 of the electroproduction process 

e(h) + N( Pl ) e(fc 2 ) + ttV) + V (p v ) + N(p 2 ) (39) 

or 

7 *(«) + N( Pl ) -> 7r°( P7r ) + niPv) + N( P 2). (40) 

To perform a leading order computation of such process (see Fig. 0}. we need to introduce the concept 
of generalized distribution amplitude (GDA) 31 for 71-77. 



4.1 7Y7] generalized distribution amplitude 

In this subsection, we briefly introduce and discuss the generalized distribution amplitude related to the 
7rfy-to-vacuum matrix element. On the basis of Lorentz invariance, the Tr a i] GDA may be defined 9 as : 

^ (^)^b i; )|^/ 2 ("^/2)7 AI h^/2;z/2]T / 3 2/i ^ /l (-z)|0) = 
1 

p£„ J dye^-y^/^n) ^ C m *j % (41) 
n 

where the total momentum of irr] pair is = p^+Prj while m 2 .^ = p 2 v . We omit the Q 2 dependence of 
the ir°r] GDA's which is the same as the one discussed above for the hybrid meson distribution amplitude. 
Note that the ir-q distribution amplitude <i>( 7rl ') describes non resonant as well as resonant contributions. 
It does not possess any symmetry properties concerning the ^-parameter. 

Let us now discuss the £ parameter. When the two mesons have equal masses, the parameter £ is usually 
defined as C = pt /p + - In the case of two different particles it is more convenient to define the parameter 

S A very similar analysis may be carried for the nr]' decay mode of the candidate 7ri(1600). 
9 A straightforward generalization enables to write similar equations for charged states 
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Figure 4: Typical diagram describing the electro-production of irrj pair. The higher and lower blobs 
represent the GDA's and GPD's, respectively. 

( in the following way: 

z _ pt ml -ml 



{Ptt+P v ) + 2ml v 

A- 2 9 

W= (p*+P„) + + 2ml v ■ (42) 

Then, we get the ordinary relation between £ and the angle 9 cm , defined as the polar angle of the tt meson 
in the center of mass frame of the meson pair: 

2C- 1 =/3cosfl cm . (43) 

In l|43(l , the standard 0- function is given by 

/3 = ^L, (44) 
m-xri 

where |p| denote the modulus of three-dimension momentum of tt and 77 mesons in the center-of-mass 
system. 

In the reaction under study, the irrj state may have total momentum, parity and charge-conjugation in 
the following sequence 

J PC = 0++ 1-+ 2++, ... 
that corresponds to the following values of the nr/ orbital angular momentum L: 

L = 0, 1, 2, ... 

respectively. We can see that a resonance with a irr/ decay mode for odd orbital angular momentum L 
should be considered as an exotic meson. 

The mass region around 1400 MeV is dominated by the strong a2(1329) (2 ++ ) resonance It is 

therefore natural to look for the interference of the amplitudes of hybrid and ai production, which is 
linear, rather than quadratic in the hybrid electroproduction amplitude. Such interference arises from the 
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usual representation of the 71-77 generalized distribution amplitude in the form suggested by its asymptotic 
expression : 



$^)'°(y,C;m^) = 10y(l-y)cf /2) (2y-l)^B li (m^)P i (cos 



(45) 



Keeping only L = 1 and L — 2 terms, we model the ttij distribution amplitude in the following form: 



$^)(y,C,m^)=30y(l-y)(2j/-l) 



Bniml^Pi (cos 6) + B 12 (m^)P 2 (cos ( 



(46) 



with the coefficient functions i3n(mj r/ ) and B^im 2 .^) related to corresponding Breit-Wigner amplitudes 
when m 2 .^ is in the vicinity of M 2 2 , M^. We have (see the technical details how to calculate these 
coefficient functions in Appendix A): 



Bii(m^) 



, 2 „ M ,~ 3 M\ -ml v -iT H M H 

7TT7 H 



and 



10 i9a 2 ^fa 2 Ml(3 2 
9 Ml-ml n -iT a2 M a2 - 



(47) 



(48) 



In the a-i case, we use the results and conventions of [33]; note that the coupling constant g a2 nri nas mass 
dimension equal to —1. 

The coupling constants ga^r) and g a2 itr) may be estimated through the approximate measurements of the 
partial widths of the a 2 and hybrid meson in the Ttr) decay channel, we have: 



r(a 2 - Trr?) = ^gl 2 „ 



\*{M 2 H ,ml,ml) 



M 5 H 



X 5 (M 2 2 ,ml,m 2 ) 



M 7 

a 2 



Neglecting the masses of n and rj mesons compared to the a 2 and hybrid meson masses, we get 



2 _ 167r rvw \ 2 
9Hw n W ^ nV ^' 9a27rv 



24tt 

lr± a 2 



r(a 2 -> ttt?). 



(49) 



(50) 



4.2 Differential cross section for in] electroproduction 

Let us first fix the kinematics. For reaction (|39|l we choose: 

fci = (ei, £1 sin 0i, 0, eicos^i) 

Pl = {E u 0, 0, p\), q = h ~ k 2 = (g , 0, 0, -pf), 
p 2 = {E 2 , |p»2 1 cos</>sin#, |p 2 | sin^>sin0, |p 2 |cos#). 



(51) 



The following Mandelstam and dimensionless variables can be defined as (here, the "hatted" symbols 
refer to the subprocess (001) 



s = {prvr, +P2) 2 = {q + Pif, i=(p 2 -pi) 2 , S=(ki+piY 



(52) 
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N(p2) 



Oz* 




N(pl) 



Figure 5: Typical process describing the electroproduction of a irr] pair 



and 



xb 



pi ■ q 



2p\-q pi -ki 

In l|51|) . the energies and momenta can be expressed as 10 



£l = 
Ei 



S - m 2 N - Q 



2 VI 



£2 = 



s + m 2 N + Q 2 



2V§ 



2VI 

E - ^ 



s ~ ml„ + m 2 N 



s — m 



Pi 



X(s,m%, 



2 ) 



2a/I 



2a/I 

Pt^I = |P2| 



9o 



2 - O 2 



2Vs 
\{s,ml v ,m 2 N ) 



2Vs 



where the kinematical function A in defined in (|30|) . 
The corresponding angles take the forms 11 

2Q 2 s 



s — m 



cos tf 2 



cos 9 — 



N Q 



(S-s)X(s,m 2 N ,-Q 2 ) X(s,m 2 Nl -Q 2 Y 
2s{t-2m 2 N ) + (s + m 2 N + Q 2 )(s - m 2 + m 2 N ) 



X(s, m 2 N , -Q 2 )X(s, m\ v m 2 N ) 
It is useful to note the following relations between the invariants: 



xb = 



Q 2 



s + Q 2 - m 2 N ' 



Vl 



Q 2 



s + Q 2 - m 2 N 



x B {S-m 2 N ) S-m 2 N 



x B yi(S - m 2 N ), s 



10 Here, 62 is the energy of the scattered lepton. 
xl (92 defines the polar angle of the final lepton. 



1 - x B 
xb 



Q 2 + m 



N ■ 
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One may also work within the center-of-mass system of the meson pair, where we have after the corre- 
sponding boost, 

Ptt = (E„, |p|sin0 cro , 0, |p|cos6> cm ), p v = (E v , -|p|sin0 cm , 0, -|p|cos0 cm ), (57) 
where the energies and momenta of the mesons take the forms 



E„ 



ml„ - ml +ml _ ml„ - ml + ml , \(m\ m 2 m\ ) 



■KT) rj 1 7T 7T7/ 7T 1 7/ 



2m T: 



2m T: 



2m„ 



(58) 



We now come to the expression for the differential cross section of reaction l|39|) . The amplitude of this 
reaction is given by 



T 71 " = u(k 2 , s 2 )i ■ e L u{ki, si)-^Aj q) v , 



leading, keeping the leading order in Q 2 , to 



T (<z) 



The amplitude of subprocess l|40(l reads 

A 



' 1 _ Q 2 yf IA 
w ° n eTTa s CF 



(59) 



(60) 



(9) " N C Q 

B 11 (ml n )P 1 (cos9 cm ) + B 12 (m^)P 2 (cos 9 cm ) 

Finally, the differential cross section of process (|39f) takes the form 

1 



da n<>r i 



dQ 2 dyididm. Kri d(cos9 cm ) 4(4tt) 5 yi\ 2 (s,-Q 2 ,m 2 N ) 



rpiT 77|2 



(61) 



(62) 



5 Calculation of the angular asymmetry 



Asymmetries are often a good way to get a measurable signal for a small amplitude, by taking profit 
of its interference with a larger one. In our case, since the hybrid production amplitude may be rather 
small with respect to a continuous background, we propose to use the supposedly large amplitude for a 2 
electroproduction as a magnifying lens to unravel the presence of the exotic hybrid meson. Since these 
two amplitudes describe different orbital angular momentum of the tt and rj mesons, the asymmetry which 
is sensitive to their interference is an angular asymmetry defined by 



/ c os 9 cm da^ " 1 (Q 2 ,yi,i, m,, , cos 9 c m ) 
J da 7T ° r i(Q 2 ,yi,t, m^cos 9 cm ) 



A[Q ,yi,t,m vn ) = ; 7. — n — (63) 



as a weighted integral over polar angle 9 cm of the relative momentum of it and r\ mesons. The angle 
9 cm is related to the parameter £ by formula jj^l . Due to the fact that the cos # cm -independent factors 
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Angular Asymmetry 
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AA 
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a' 
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k 



Figure 6: The angular asymmetry as a function of m 



in both the numerator and denominator of l|63|) are completely factorized and, on the other hand, these 
factors are the same, we are able to rewrite the asymmetry as 



A(m nv ) 



J d(cos9 CHI ) cos( 



En (ml )Pt (cos 6 cm ) + Bi 2 (ml )P 2 (co$6 cm ) 



J d(cos 9 C 



B 1 i(m2 ?) )P 1 (cos6' cm ) + B 12 (ml n )P 2 (cos 6 cm ) 



7T7] ) 



(64) 



and, calculating the cos # cm -integral analytically, to obtain 



A(m Tr ,) 



N(m VTi 



7T7] ) 



with 



B ii(ml v )Bl 2 (ml v ) 



Bi 2 (m 



2 > 

TTTjy 



(65) 



(66) 



While in the two-pion production case the interference between the isoscalar and isovector channels can 
be investigated, we here restrict to the interference between L = 1 and L = 2 modes of 7T°r]. As a result, 
the introduction of the so-called intensity density (see [H]), i.e. the integrated-over- invariants value, 
useful for the two-pion modes, completely coincides with the value l|65(l for our case. 
Our estimation of the asymmetry l|65|l is shown on Fig. Since the numerator of 1)65(1 , i. e. the real part 
of the product of B\\(m\^) and B* 2 (ml v ), is proportional to the cosine of the phase difference ASi^ — 
Si = i — Si =2 the zeroth value of H65|) takes place at A6i t2 — tt/2. This is achieved for w 1.3 GeV. 
Besides, one can see from Fig. [5] that the first positive extremum is located at m^^ around the mass of 
a 2 meson while the second negative extremum corresponds to the hybrid meson mass. 
Note that this angular asymmetry is completely similar to the charge asymmetry which was studied in 
electroproduction at HERMES 35 . 
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6 Note on the 7T7r7r channel 



The hybrid candidate 7Ti(1600) has been also seen through a nun decay channel. In that case the deep 
exclusive electroproduction of three pions provides the background which should be studied, including the 
possible interference effects of hybrid meson and the background. The analysis that we have described 
in section 4 and 5 may be adapted to the three body case by using the results of At the leading 
twist level the generalized distribution amplitude of the charge conjugation even state may be written in 
complete analogy to the pion light-cone distribution from the large distance matrix element 

SaP = ^J dx ~ e ~ iv[P+x ~ ] ou>H M^-v')Mo) I0) in (67) 

as 

S 9>a p 7^75 = f Co, C+, C-5 W 2 12 , , W* , ) P+. (68) 

The three light cone fractions are normalized by the condition Co + C+ + C- = 1, making only two of them 
independent, while the squared total energy of the three pions is W 2 — W 2 2 + W 2 3 + — 3mJ, where 
Wij are the invariant masses of the pairs of mesons. 
Charge conjugation invariance provides a symmetry relation: 

$ + (y, Co, C+, C-) = - V, Co, C-, C+). (69) 
The asymptotic z-dependence is just 

$+(!/, Co, C+, C-) = ^TajVi 1 ~ V)«0 + aCo 2 ), (70) 

where a is an unknown parameter and the normalization has been hxed with the help of the fact that 
putting both charged pion momenta to zero, one should get the GDA equal to the pion distribution 
amplitude. The QCD evolution is the same as for the pion distribution amplitude, i.e. with a vanishing 
anomalous dimension. 

Conversely, the generalized distribution amplitude $ _ (y, Co, C+, C-) °f the charge conjugation odd state 
obeys the equation : 

$- (y, Co, C+, C-) = - y, Co, C-, C+). (7i) 

and its asymptotic z-dependence is 

$-(y, Co, C+, C-) ~ 2/(1 - »)(2y - l)P(Co, C+, C-), (72) 
where P is a polynomial of degree 3 with the symmetry property 

P(Co,C + ,C-) = J P(Co,C-,C+)- 

Its QCD evolution is the same as for the hybrid distribution amplitude. The interference with 7Ti(1600) 
may produce an angular asymmetry similar to that of Section 5, although the appearance of third pion 
makes the analysis more complicated. 



7 Diffractive production of the hybrid at large energy 



At large energy, one should consider a different framework, namely the impact representation 2 where 
the meson electroproduction amplitude is factorized in impact factors and a reggeized two or three 
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gluon exchange known as perturbative Pomeron or Odderon exchange. Without performing a detailed 
phenomenology of this reaction in this regime let us recall some well known formula and briefly propose 
a strategy to help future (or present) experiments at high energy to search for hybrids. 
The even charge conjugation of the hybrid meson selects in this case the Odderon exchange [37] and the 
amplitude for its electroproduction is equal to 



Mn = — 



8tt 2 s f d 2 kid 2 k 2 d 2 k3$ (2) (ki + k2 + k 3 -p ff o) 



3! 



(27r)6k?k|k§ 



77 



J, 



O 



(73) 



where Jq ^ H (ki,k 2 ,k 3 ) and Jq~* n (ki,k2,k2) are the impact factors for the transition 7* — > H° via 
Odderon exchange and of the nucleon in initial state N into the nucleon in the final state N'. The 
impact factors are defined as the s— channel discontinuities of the corresponding S— matrices describing 
the 7*0 — * H° and NO — > N' processes projected on the longitudinal (nonsense) polarizations of the 
virtual gluons in the t— channel. 

The upper impact factors are calculated by the use of standard methods. The leading order calculation 
in pQCD gives in the case of a longitudinal polarized photon : 



j£{\i 1 MM) = - ie \ d ^ Q J d V yyPo^xMM)\^ H °{y) 



(74) 



where ki + k 2 + k 3 = p ff o and 



fb(ki,k 2 ,k 3 



9 ? 

JTPjjo 



1 



3 



V 2 Pho 



(kj - yp H o ) 



- /'- 



(75) 



The proton impact factor cannot be calculated within perturbation theory. One may use phenomenolog- 
ical eikonal models of these impact factors proposed in Ref [35| which read 



TiV-tiV' -9 a 



3 Jabc 



4iV c 



3 

F(p H o , 0, 0) - ^ Ph" - ki, 0) + 2 F(ki, k 2 , ka) 

i=l 



(76) 



where 



and A 



^(ki^ka) 



.4 2 



+ i [(k a - k 2 ) 2 + (k 2 - k 3 ) 2 + (kg - ka) 2 ] 



(77) 



In these equations we denote the soft QCD-coupling constant by g and one may take 
o<-soft — 3 2 /(47r) = 0.5 as a reasonable mean value. 

Since the Odderon amplitude is known to be rather small, producing the hybrid in electroproduction at 
large energy will be rather difficult. It will thus be useful to search for the hybrid meson in this context 
through an interference of a Pomeron mediated amplitude to a Odderon mediated one, as discussed in 
PU). The three pion channel discussed in the preceeding section is interesting in this respect since a 
charge asymmetry between the 7r + and the 7r _ will single out this interference. The charge asymmetry 
to measure is defined as a integral weighted with an antisymmetric function in the exchange C, + 
the simplest example being ( + — £ _ . : 



A(Q 2 ,i,m 2 J 



f d( + d(- {(+-(-) 2 Re Ml L (M L y 
L 







\Mf\ 2 



'o 



(78) 
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which is approximately equal to the ratio of the Odderon exchange amplitude and the Pomeron exchange 
amplitude, which one may approximate to the production of a J PC = 1 state like w(1650). 
This may be related to a forward-backward asymmetry in the rest frame of the two charged pions, the 
forward direction being defined as the direction of the neutral tt° meson. Defining the angle as the 
angle of the ir + to the 7r° 3-momenta in this frame, one gets 



A(Q 2 ,t, 



J cos 9 da(s,Q 2 ,t,'i 



/ da(s,Q 2 ,t,ml„, 



2 / d(cos6) cos 6*3? M£(Mg) 
-l L 

1 r * » _ 

/ dcosO \Ml L \ 2 + \M 7 L \ 2 



(79) 



8 Conclusion 



In conclusion, we have calculated in this paper the leading twist contribution to exotic hybrid meson with 
J PC = 1 ^ electroproduction amplitude in the deep exclusive region. The resulting order of magnitude 
is somewhat smaller than the p electroproduction but similar to the tt electroproduction. The obtained 
cross section is sizeable and should be measurable at dedicated experiments at JLab, Hermes or Compass. 
We made a systematic comparison with the non-exotic vector meson production. To take into account 
NLO corrections, the differential cross-sections for these processes have been computed using the BLM 
prescription for the renormalization scale. In the case of p production, our estimate is not far from a 
previous one which took into account kinematical higher twist corrections. 

We have also discussed in detail the ttt] mode corresponding to the 7Ti(1400) candidate in the reaction 
ep — > epir a r]. We have calculated an angular asymmetry implied by charge conjugation properties and 
got a sizeable hybrid effect which may be experimentally checked. 

In the region of small Q 2 higher twist contributions should be carefully studied and included. Note that 
they have already been considered in the case of deeply virtual Compton scattering 0D] where their 
presence was dictated by gauge invariance, and for transversely polarized vector mesons |41| where the 
leading twist component vanishes. We leave this study for future works. 

Finally, the diffractive production at very high energy has been briefly studied. The weakness of Odderon 
mediated processes makes the study of hybrid meson electroproduction a very difficult task for HERA 
experiments. 
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Appendix A: Functions Bn(mD and B^m^) 



TTT/J 



We now proceed to the calculation of the functions Bu(mL) and i? 12 (m^. j; ) related to the corresponding 
Breit-Wigner amplitudes. Let us start from the consideration of the Trry-to-vacuum matrix element of 
some vector nonlocal quark operator. We have 



(7T(p^ V (Pn)\Ol(~z;z)\0), 
where = ?/>(— 0/2)7^(2/2). This matrix element can be rewritten in the equivalent form: 



(Tr(p„)r,(p v )\H(p))- 



1 



(80) 



(81) 



u* ■ r -Tl-(H(p)\O^(-z;z)\0) + (other reson.). 
Mfj-p 2 -iY H M H 

We will, from now on, neglect the contribution from other resonances. Note that owing to the momentum 
conservation law we have p 2 = (p^ + p n ) 2 = w„„. 

Further, we introduce the parameterization of the relevant matrix elements: 



(n(Pir)'n(Pri)\H(p)) = G { H l 71 {ml,m 2 v , Mjj)(p^ - p v ) ■ e (A) = -igu-Kr,^ ~Pv) ' e 

(H{p)\0^(-z;z)\Q)=if H M H - p^ / dye^~ 2 ^l 2 <^(y\ 

p ■ z j 



(82) 



Using the above-mentioned parameterization, equations QSOjl and (|81|l may be rewritten as 



(P.+P V )J dye^-^+^'^^iyZrn 2 ^) = 


i/hMh {-i)9H-wnPp (Ptt - Pvf 772 — 7^2 r U 

p-z Mfj - mi v - iT H M H 



dye 



i(l-2y)p-z/2,H 



(83) 



The summation over polarization vectors reads 

/ , e a C /3 — 



therefore the contraction of corresponding vectors with (|84l) gives us 



(Ptt -P n ) + , m l~m 2 



Ml 



(84) 



(85) 



p • ^ ■ ■ v y (ptt +^) + 

Further, the term B\\ corresponding to the hybrid meson, see l|46|l . can be rewritten in the form: 
^4(y,C,ml v ) = 18y(l-y)(2y-l)B 11 (m^)P 1 (cos^), 

where eqn. i|43[l is used. Inserting this function into the Ihs of l|83|) we get an explicit expression for the 
function B\\\ 



(86) 



18 5 Mfj - m 2 n - iT H M H 



_ mj - m 2 ( Pn + Pr] )+ 



mi-m 2 / ( Pn - Pn )+ to; - m; n 



M 2 H 



(87) 
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where 



T 



dye * (1 - 2yW 2 H (y)) 



If we use the asymptotic form for the function (p^ (y), the expression (|87|l is rewritten as 



3 Ml 



ml n - iT H M H 



I - m l ( (Ptt - Prj) + 

(Pn +P n ) + 



71 Tj 



(89) 



where the value m 2 ,^ is in the vicinity of the hybrid mass Mh- We thus have obtained the following 
expression for the function Bu 



Ml 



3 Ml 



mlr, - iT H M H 



hi-, 



(90) 



Ml 



We will now focus on the calculation of the function Byiijn^^. As mentioned above, the a2-resonance 



formula for the case 



Ml 2 reads (cf. P|l and ijST}) 

(7r(p„)v(p n )K(-z;z)\0) = 

(*(P«) v(Pv)\02(p)) M2 _ 2 V r M (a 2 (p)\OX(-z; z)\0). 



(91) 



The parameterizations of matrix elements standing in (|91|) can be introduced in the following forms. First 
of all, we write the parameterization for the in]-to-a2 matrix element: 



W^)r ? (p I7 )|a 2 (p)) = eWv^, 



where 



0, 



0. 



(92) 



(93) 



(94) 



Due to the Lorentz invariance the most general representation of the tensor V^ v take the form 

V"" = G%l v {ml,mlMl)p»pl + G«„(m*, m% M 2 2 )p^ + 
G& nv (ml , m% , M 2 a2 )p£j£ + G^ 2 \ n (mj , m* , M 2 2 )p%p» . 

We fix the form factors as 

Ga2TTTi( m TTJ m r]j M a2 ) = G^ 2 \ v (ni-n , rn^ , M a2 ) = ig a2 iTn' 

Ga 2 \r] ( m ir i m r) ; ^a 2 ) = ^a 2 Trrj( m in TO r/? ^a 2 ) = ^Wa^rj- 

In this case the parameterization of (|92|l is reduced to the form 

(""(Ptt) r](p v )\a 2 (p)) = ig a , 2 wT,e'$ (Ptv - P V Y (p* - P V Y- (96) 



(95) 
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The parameterization of vacuum-a2-meson matrix element can be written as 



(a 2 (|,)|^(-z;z)|0)=/ a2 M ( 2 ' 



P/1 dye^-^/^iy). 



{p-zf 

u 

Again, the asymptotic form of function <jy£ (y) can be defined as 

$ a L *(y) = 30y(l-y)(2y-l). 
Further, the a2-resonance part of nr] distribution amplitude reads 

$%${y, C, ml v ) = 182/(1 - y)(2y - l)S 12 (m' )P 2 (cos 



(97) 



(98) 



(99) 



Thereafter, inserting the parametrical representations of hadron matrix elements l|96() . (|97() in eqn. I|91[) 
with (|98|l . Ij99(l . one can see that 



B 12 (ml n )P 2 (cos 9) = 

30 i9a2TTyfa2 M l 2 

18 Ml - m|„ - iT a2 M a2 



(100) 



z a z 13 



(p ■ z) 2 



E 



P (A) *(A) 



where as usual 



Ve (A) e* 

*V e «/3 
A 



Pa\Pa 2 

M 2 

a 2 



Using 143(1 , the Legendre polynomial P 2 (cos 9) can be presented as 



P 2 (cos( 



- cos t) = 

2 2 



(101) 



(102) 



7TT) 



2\ 2 (ml r) ,m 2 ,m 2 ) 
A straightforward computation leads to 



3(2C - l) 2 - 



\ 2 {ml n ,m 2 ,ml) 



711) 



^ z a z? ^ j X) j {x) _ 1 2P 2 (cos9)X 2 (ml J1 ,m 2 ,ml) 

mi 



3 (3(2C-1) 2 



where we used that 



(jp.z) 
X 2 (ml v ,m 2 ,ml' 



2 n M l 2 -ml-m 2 
z =0, p 7I -p, 1 = . 



7H/ 



(103) 



(104) 



We have finally obtained the following representation for the function J5i 2 (m 2 r?7 ) in the vicinity of m 2 .^ 

iga 2 7r V fa 2 M 2 2 X 2 (ml n ,m 2 ,ml) 



M 2 ■ 

a 2 



10 

Ym 2 , 



a 2 ~ ' n lr, ~ iT o, 2 M a2 



(105) 
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